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Consider  the  equation 


(1)  ut  =  (aux)x  +  f(u),  0  <  x  <  1, 

with  homogeneous  Neumann  conditions 

(2)  =  0  at  x  =  0,  x  =  1, 

2  2 
where  a  £  C  [0,1]  is  a  positive  function  and  f  £  C  QR) .  System  (1),  (2) 

1  12 

defines  a  local  dynamical  system  in  H  (0,1)  =  W  *  (0,1)  (see  Henry  [8]). 
Furthermore,  the  co-limit  set  of  any  bounded  orbit  is  exactly  one  equilibrium 
point;  that  is,  a  solution  of  the  equations 

(3)  (aux-*x  +  f^U')  =  °’  0  <  x  <  1 

(4)  u  =  0  at  x  =  0 ,  x  =  1 

x 

(see  Matano  [11],  Hale  and  Massatt  [6],  Zelenyak  [15]). 

The  first  result  is  the  following. 

Theorem  1.  if  aM(x)  £  0  on  [0,1],  then  every  nonconstant  equilibrium 
solution  of  (1),(2)  is  unstable. 

For  a(x)  =  c  >  0,  where  c  is  a  constant,  this  result  was  proved  by 
Chafee  [2].  Note  that  the  conclusion  in  the  theorem  is  valid  for  every  function 
f.  Matano  [13]  has  given  an  example  (1),(2)  for  which  there  are  stable  non¬ 
constant  equilibrium  solutions.  The  function  f  can  be  chosen  to  be  a  cubic 
with  three  simple  zeros  and  the  function  a  is  1  on  intervals  [0,ot],  [6,1] 
and  <_  e  on  [y,<5],  a  <  y  <  6  <  ft  and  t  is  sufficiently  small.  Numerical 
examples  indicated  the  same  property  can  occur  with  a  >_  1  on  [0,a]  and  <_  e  on 
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[8,1]  with  a  <  8.  These  examples  show  clearly  that  more  work  is  needed  to 
determine  the  optimal  conditions  on  a  for  the  conclusions  of  the  theorem  to 
hold. 

For  a  more  general  parabolic  equation  in  several  space  dimensions  in  a 
bounded  domain  ft  and  the  diffusion  coefficients  constant,  Casten  and  Holland 
[1],  Matano  [12]  have  shown  that  Theorem  1  is  true  if  ft  is  convex.  The 
authors  have  had  no  success  in  extending  this  result  in  a  significant  way  with 
variable  diffusion  coefficients.  If  ft  is  not  convex  and  the  diffusion  co¬ 
efficients  are  constants,  nonconstant  equilibrium  solutions  can  occur  (see 
Matano  [12],  Hale  and  Vegas  [7],  Vegas  [14])  for  certain  nonlinear  functions  f. 

If  the  diffusion  coefficient  in  (1)  is  constant  and  f  is  a  function  of  u(-,t), 

r  1 

for  example,  f(u(x,t))  =  g(u(x,t),  a(x)u(x,t)dx) ,  then  stable  nonequilibrium 

J0 

solutions  can  occur  (see  Chafee  [3]).  A  similar  situation  was  considered  in 
several  space  dimensions  and  nonconvex  domains  by  Keyfitz  and  Kuiper  [9].  For 
constant  diffusion  in  (1)  and  f  replaced  by  s(x)f(u),  Fleming  [5],  (see  also 
Henry  [8])  has  considered  how  the  existence  of  stable  nonconstant  solutions 
depend  on  s(x).  All  of  these  papers  should  be  reconsidered  with  variable 
diffusion. 

There  is  an  analogue  of  Theorem  1  for  the  Dirichlet  problem. 

Theorem  2.  Consider  the  equation  (1)  with  Dirichlet  boundary  conditions 

u  =  0  at  x  =  0,  x  =  1. 

If  a"(x)  ±  0  on  [0,1]  and  v(x)  is  a  nonconstant  equilibrium  solution  such 
that  vx  =  0  at  two  points  in  (0,1),  then  v  is  unstable. 
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For  a(x)  =  c  >  0,  this  result  was  proved  by  Chafee  and  Infante  [4], 
Maginu  [10]. 


Proof  of  Theorem  1 .  If  v  is  a  solution  of  (3) , (4)  and 

rl 


= 


[a *1  -  f'(v)<t>2]dx, 

\  A 


then  the  first  eigenvalue  X^  of  the  operator 


Lu  =  -  f '  (v)u 


is  given  by 


X  =  :  *  e  H1(0,1),|  +  |  ?  =1}. 

L  (0,1) 

Furthermore,  the  equilibrium  solution  v  of  (1),(2)  is  unstable  if  X^  <  0 
(see,  for  example,  Henry  [8]).  Consequently,  it  is  sufficient  to  show  that 
the  hypotheses  of  the  theorem  imply  X^  <  0  if  v  is  not  a  constant  and  this 
will  be  the  case  if  we  show  that  -<  0  if  v  is  not  a  constant.  We 

have 


#"(avx) 


£ 

£ 


[a[(avx)x]2  -  f' (v) (avx)2]dx 


af  (v)dx  - 
0  J 


1 


f(v)  a  av  dx 

X  X 


[  af2(v)dx  +  [  f 
J  0  J  0 


af  (v)dx  +  J  f (v)a(avx)^dx  + 


£ 


(avx)xa' (avx)dx 


f(v)aTav  dx 

J0  x 
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To  compute  the  last  expression,  we  integrate  by  parts  to  obtain 

rl 

av  af(av  )  dx 

X  XX 


fl 

’I 

(av  )  a’(av  )dx  =  - 

av  aMav  dx  - 

JQ  x  X  x  J 

Iq  X  x  J 

and,  thus. 


Therefore, 

C5) 


i: 


(av  )  a1 (av  )dx  =  - 

XXX 


2 

a"(av  )  dx. 
0  X 


^'(avx)  =  ± 


a"(av  )  dx  £  0 

A  X 


since  we  have  assumed  a11  £  0. 

If  _£#(av^)  <  0,  then  X^  <  0  and  the  solution  v  is  unstable.  Thus, 

suppose  X^  >_  0.  From  (5),  this  implies  =  0  which  implies  X^  =  0. 

Since  the  first  eigenvalue  of  L  is  simple,  there  is  a  <j>  satisfying 

L(j>  =  0,  <J>  =0  at  x  =  0,  x  =  1,  |<f>|  9  =1.  Also,  ^(av  )  =  0  implies 

X  L  (0,1)  X 

there  is  a  constant  c  such  that  av  =  c<f>  on  [0,1].  If  c  ^  0,  then  v  =  0 

x  x 

at  x  =  0,  x  »  1  and  a  >  0  imply  <(>  =  0*  at  x  =  0,  x  =  1.  But  this  would 
imply  4>  =  0  on  [0,1]  which  is  a  contradiction.  Thus,  c  =  0  and  av^  =  0, 
v^  “  0  and  v  =  constant.  This  proves  the  theorem. 


Proof  of  Theorem  2.  Suppose  v  is  an  equilibrium  solution,  vx  =  0  at  x  =  a, 
x  =  3,  a, 6  €  (0,1).  Then  v  is  an  equilibrium  solution  of  the  Neumann  problem 
on  the  interval  [a,  3].  As  in  the  proof  of  Theorem  1,  one  shows  that  the  first 
eigenvalue  X^  of  the  linear  variational  operator  is  negative  if  v  is  not  a 
constant.  By  the  characterization  of  X^  as  a  minimum  of  the  functional  &(u) , 
it  follows  that  X^  for  the  Dirichlet  problem  on  [0,1]  is  negative.  This  proves 


the  theorem. 
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